Chern numbers and localization by non-Hermitean operators 
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In Hermitean quantum mechanics, extended current-carrying states are distinguished from localized 
ones by their non-zero Chern number. We generalize the notion of Chern number to non-Hermitean 
localization problems such as tilted flux lines in type-II superconductors with line defects and passive 
scalar transport in random vorticity fields. We find that the usual point eigenvalue degeneracies that 
occur in the Ifermitean case are now loops bounding a branch cut sheet with the same quantized 
total Berry flux, and hence Chern number, as the original point. If a loop cuts the integration 
surface non-quantized contributions to the total flux result. 
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In Hermitean quantum mechanics the differential geo- 
metric notion of Chern number arises in the context of 
adiabatic transport of eigenstates in Hilbert space 0. 
The change in the overall phase of a state after it is 
transported around a closed loop contains a purely ge- 
ometric contribution associated with curvature of the 
Hilbert space. Quantized integer Chern numbers are as- 
sociated with singularities in the curvature that originate 
in eigenvalue degeneracies of the Hamiltonian generating 
the transport. 

One way of transporting an eigenstate is by changing 
its boundary conditions. A Chern number can then be 
computed as a certain integral over the space of bound- 
ary conditions. Since a localized state is insensitive to 
boundary conditions, a nonzero Chern number should be 
associated only with an extended state, and serves as a 
nontrivial diagnostic of localization. This procedure has 
been used to study localization in the integer quantum 
Hall effect ||^ . Physically, the Chern number reflects the 
ability of the state to transport charge across the system, 
and is proportional to its Hall conductance 

Recently there has been much interest in localization 
properties of certain non-Hermitean operators, for exam- 
ple, operators governing effective two-dimensional quan- 
tum descriptions of pinning of tilted flux lines by ex- 
tended defects in three-dimensional type-II superconduc- 
tors , and passive scalar transport by two-dimensional 
random vorticity fields Q. Both these examples have 
close ties to the problem of electron localization in ran- 
dom magnetic fields ||^. In this paper we will consider 
the generalization of Chern numbers to non-Hermitean 
operators. We review the non-Hermitean Berry phase |^] 
and how it gives rise to Chern numbers through eigen- 
value degeneracies. Application is then made to a class of 
operators that interpolate, as a function of a parameter 
< 6 < tt/2, between a special case of the Hermitean 
random flux model and the non-Hermitean flux line and 
passive scalar problems. 

Consider a family L(X) of (possibly non-Hermitean) 
operators, where X = {Xi,X2, ■ ■ ■) is a set of fixed real 



parameters spanning a manifold. Let |m;X) and (m;X| 
be the right and left eigenvectors of L(X) satisfying 



L(X)|m;X) 
(m;X|L(X) 
(to; X|n; X) 



Era{X){m;X\ 



(1) 



The eigenstates are determined up to an arbitrary com- 
plex amplitude, and we make a fixed arbitrary choice 
of this quantity so that |to;X) is a single-valued func- 
tion of X over the manifold. The projection operator, 
P„j(X) — |r7i; X) (m; X|, is independent of this ampli- 
tude. The eigenvalues i?rn(X), assumed nondegenerate 
for now, may be complex, and in general (TO;X|t y^ 
|m;X). For any given operator L(i), we define the dy- 
namics of states |(/)(t)) and {(j){t)\ via 

tdtm)) = mm)), - ^dtim = (mm- m 

This dynamics preserves all inner products (0i(i)|02(t)). 

Let X(a), < a < 1, X(0) = X(l), parameterize a 
closed path, C, in parameter space. We define parallel 
transported eigenstates \m;a) = 6m(a)|TO; X(a)) along 
the path X(ck) by the rule 

dabyn{a) = -{m;X(a)\da\m\yi{a))bm{a), (3) 
with 6,„(0) = 1. The solution is ^ 

.X(a) 



i\ A„,(X).dX 

C:X(0) 



hm.{a) = exp 
A„(X) = (TO;X(a)|i(9x|TO;X(Q:)) 



(4) 



Direct substitution shows that this parallel transport 
may be generated dynamically using the operator L(<) = 
i[5tP„(X(t)),P™(X(i))] (with a = t) §. For a closed 
loop one obtains 



&™(1) = e*^"^''\ 7m(C) = / dX • A™(X). 

Jc 



(5) 
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In the Hermitean case A™, and hence 7(C), is real. In 
the non-Hermitean case they are both generally complex. 
It is easy to see that 7(C) is independent of the arbi- 
trary normalization of the states |m;X). Any change 
of normalization, e'^^"^-* , appears as a gauge transforma- 
tion A' = A — 9xA, under which the integral in (||) is 
invariant. Provided no singularities are encountered we 
may write 7(C) as the integral of the gauge invariant flux 
using Stokes theorem, 

F™(X) = V X A„(X) 

= ^[ax(™o;X|] X [ax|mo;X)], (6) 



over any surface S bounded by C 



lm{C) = / dI]-F„(X). 

Js 



(7) 



By construction one has V • = 0, as is required for 
7(C) to be independent of the choice of S. 

By inserting a complete set of states, the flux may be 
written in the form 



F„,„(X)=z5] 



,(m;X|9xL|n;X) x (n; XI^xLIto; X) 



[En — 



(8) 



where the prime indicates that the term n = m is omit- 
ted. Equation ^ demonstrates explicitly that so long as 
L(X) is smooth, singularities in Fm can only occur when 
two eigenvalues cross: En — Em — s- for some n ^ m. 
In the Hermitean case eigenvalue crossings, and there- 
fore singularities in Fm(X), generically occur at isolated 
points in the three dimensional space of X. One may 
choose the surface S in (|^) to avoid such points by pass- 
ing either over or under them. On the one hand, the 
difference of the integrals over two such surfaces 5*1 and 
S2 is 



lm{S) = 7m(S'l) 



Js 



(9) 



where 5 = 5i U 6*2 is a closed surface enclosing the singu- 
larity. On the other hand, e'^™^'^^^ = e'^™^'^^^ since both 
surfaces are bounded by the same curve C. We conclude 
that for any closed surface 5, 7™ (5') = 2T:q for some in- 
teger Chern number q ^ . Eigenvalue degeneracies act as 
a quantized point sources of flux for F^. 

This argument does not rely directly on Hermiticity, 
but only on the result that energy level degeneracies oc- 
cur generically at points in three dimensional space, so 
that these singularities can be enclosed by singularity- 
free closed surfaces. We now address the issue of how to 
generalize this result to the non-Hermitean case. In order 
to study a degeneracy in more detail, we project L onto 
the (complex) two-dimensional subspace spanned by the 
eigenvectors |m;X) and |fT.;X). If P„j„ ~ \rn) {m\ + \n) {n\ 
is the corresponding orthogonal projection, we consider 
the matrix 



,(X) =P™„(X)L(X)P™„(X) 

= Q:o(X)I -I- Q?(X) • (T 



(10) 



where 3 are the Pauli matrices, and ao and a = 
(ai, q;2, 013) are four (at this stage arbitrary) complex 
numbers. In the Hermitean case these numbers are all 
real. The eigenvalues of L„i„ are E± = ag i AiJ, where 
Aii' = \J OL ■ OL. The left and right eigenvectors are easily 
computed as well, and the the fluxes are 



F±(X)=T 



aia2 x as -I- a2a3 x ai + aaai x a2 



(11) 



where ai(X) = 5xai(X 



Equation (11) can be used to explore the nature of the 
singularities in Fj- in the neighborhood of degeneracies, 
0. The condition for A£^ = is simply -I- -I- 
03 = 0. In the Hermitean case where the ai are all real, 
this condition requires that ol\ ^ ct2 = oiz = 0, but in 
general this equality places only two conditions on the 
six parameters (the real and imaginary parts of each a^, 
« = 1, 2, 3). In the parameter space X degeneracies will 
then generically occur on a submanifold of codimension 
2, i.e., on lines in the three dimensional space. 

Since the matrix Lm„ is no longer Hermitean, there is 
no guarantee that there will be two independent eigen- 
vectors associated with the two degenerate eigenval- 
ues. Generically Lmn will be upper-triangular (with 
a\ -\- ia2 = at the same time that as = 0, but 
ai — ia2 7^ 0), and only one eigenvector will exist, with 
the orthogonal direction spanned by a generalized eigen- 
vector |10|| . We shall call this situation a generalized de- 
generacy. What we call a true degeneracy exists only 
when both ai ± ia2 = 0, so that limn = is propor- 
tional to the identity matrix (this differs from the Her- 
mitean case only in that ao = E may be complex), and 
places six conditions on the occurence of a true degener- 
acy, which will generically not occur for a three dimen- 
sional parameter space X, unless further restrictions are 
placed on the class of matrices (such as self-adjointness). 
In what follows we do not impose any such further re- 
strictions. 



If Xg be a point of generalized degeneracy, set x 



X 



Xg. From 



11^) we obtain for small x: 



25/2(wx)3/2' 



(12) 



with complex vectors 



V = aia2 X a3 + Q;2a3 x ai + asai x a2 
w = aiai + a2a2 + a3a3, (13) 







with ai + a2 + as ~ 0. The easily proven result v • w 
guarantees that V • F± = 0. 

The condition w • x — is satisfied locally if x is or- 
thogonal to both Wfj = Rew and w/ = Imw, showing 
that the line of generalized degeneracies is locally tan- 
gent to u = w/f X w/. Note that unless it is real, v will 
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not be parallel to u. This line can pierce the integra- 
tion surface S at singular points where F± will have the 
asymptotic behavior (|l2|). Since the surface integral is 
two dimensional, the singularity will be integrable. 

Suppose now that x circles a singular line in the plane 
formed by and w/. For example, let w-x = rje^^ with 
T] fixed and — tt < 9 < tt determined by wjj ■:k — rj cos{9) 
and w/ • X = r]sin{6). Then F± cx e^*'^^/^, so that single 
valuedness can be maintained only if a surface of branch 
cuts extends out of the singular line. It is straightfor- 
ward to show that this branch sheet joins F_|_ and F_: 
tracing the closed path causes Em and En to exchange 
positions by circling each other in the complex plane. 
F± (X) forms a single complex function on a two-sheeted 
Riemann surface (see Fig. g below). 

In the Hermitean case one may deform the integration 
surface to small spheres enclosing the point singularities. 
In the generic case the deformed surface wraps around 
those portions of the lines and their branch surfaces that 
are internal to the original surface. 

In order to understand the result of the integration 
(^ we consider a model problem with the exact form 



ai • X, a 



,0 



a2 • X, a 







degeneracy exists at x = 
Near the flux is 

F±(x)=T 



when a 



■ xj, so that a true 
= (a?,aO,aO) = 0. 



2(Sx • Sx)3/2 



(14) 



and 



where S is the 3x3 matrix whose rows are the at 
Dq = det S = ai • a2 X as . 

We suppose that the a^ are all real, so that a" rep- 
resents a (generally non-Hermitean) perturbation of an 
Hermitean problem. On transformation to real coordi- 
nates y = Sx, the flux takes the classic Coulomb form 
F±(y) = TCT(£'o)y/2|y|^ where a{Do) is the sign of Dq. 
It follows immediately that the integrated flux through 
any surface S enclosing the origin is 7(5") — ^2Tra{Do). 

On including the constant shift cP, the flux becomes 



F±(y) = T^(i^o) 



a" +y 



2[(ao+y)-(ao+y)]3/2 



For very large |y| this form reduces to ([Til). The 
flux through a surface at infinity remains 7(>S'). Since 
Vy • F±(y) — 0, this result holds for any surface that en- 
closes all singularities in F±. These singularities occur at 
the zeros of (a^-l-y)^, i.e., for |y-|-Red;"p = jlma^p and 
ImcJ^ • (y -|- Rec?") = 0, which constitute a circular ring 
of radius a = |Imcf*'| formed by the intersection of the 
sphere of radius a centered at —Red" with the plane nor- 
mal to Imd" passing through this same center. In x-space 
the ring becomes an ellipse. A ring must be spanned by 
a branch surface in order to maintain single-valuedness 
of F±. 

We have therefore established that under non- 
Hermitean perturbations the usual quantized point 
charges expand into closed loops, with net flux quan- 
tized by the original Hermitean Chern number. If a loop 



grows to pierce the surface of integration, only a fraction 
of the total flux will be enclosed, and integer quantiza- 
tion is lost. This fraction depends on the geometry of 
the branch cut, which is itself arbitrary, except for its 
end points. Only the total flux, integrated over all sheets 
of the Riemann surface, remains quantized jl^. The lat- 
ter will correspond to the net Chern number of the collid- 
ing eigenvalues, which are now topologically "entangled." 
The entanglement becomes more pronounced as the non- 
Hermiticity increases, and may involve more than two 
eigenvalues that collide with each other sequentially as 
X is varied, generating further Riemann sheets. 

We illustrate this formal discussion by a speciflc exam- 
ple, a family of operators in two dimensions: 



iA{r) ■ V, 



(16) 



in which < 9 < 2tt controls the degree of non- 
Hermiticity (L is Hermitean only for 9 = 0, tt), Z? is 
a constant coefficient, and A is an incompressible vec- 
tor field: V • A — 0. For 9 — 0, this operator cor- 
responds to the random flux problem in the Coulomb 
gauge 1^, with magnetic field i? = V x A, and a partic- 
ular choice V = —A^ /AD for the scalar potential. For 
9 = t: /2 this operator represents a model of flux lines in a 
three-dimensional superconductor with extended defects, 
(where A corresponds to the horizontal components of 
the magnetic field |^), or a two-dimensional passive- 
scalar transport model, (where A is the velocity field, 
and D is the diffusion constant). We set _D = 1 and im- 
plement (16) numerically by discretizing it Q on a small 
(4 X 4) lattice. The components of the vector field A 
are initially chosen to be independent random variables 
chosen uniformly on an interval \—W, W], with W = AQ. 
Incompressibility is imposed by subtracting from A the 
function Va, where a(r) satisfies V^a = V • A with peri- 
odic boundary conditions. 

Study of random operators like (|l6|) is motivated by 
questions of localization j^-^,^. Extended states are 
distinguished from localized states by their sensitivity 
(15) to boundary conditions. Let En{4>xi4>v)j ''Pn'"'^^ {x,y) 
be the eigenvalues and cigcnstates of L on cLIl Lx -^y 
system with boundary conditions '0n^'^^(x -\- Lx^y) — 



.j<l>.^t'"'^Hx,y), i,tr*-{x,y + Ly) = e*^«^^'*'' (x, y). 
The two-dimensional space {(f>x,4'y) has the geometry 
of a torus that is embedded in a three-dimensional pa- 
rameter space X. The Berry fiux through the torus is 
the integrand of (|): Fn{4>x-,(iiy) = (^^^ Vnl^^^ Vn) - 
{dij, V'nl^^^V'n), and if the integral of the fiux over the 
surface (the Chern number) is non-zero, then the eigen- 
state is extended. 

Fig. |l] displays the variation of eigenvalues with 9 for 
(f>x — 4>y — ^- In Fig- |2| we plot the imaginary parts of 
F±{4>x,4>y) for the labeled eigenvalues in Fig. || and for 
two values of 9, revealing the emergence of a branch cut 
and showing F± are different sheets of the same analytic 
function. The log-log plot exhibits the 3/2-law ( p^ ) as 
an endpoint of a branch cut is approached. 
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eigenvector I eigenvector ll 

.e=.0253i e=.025n 



FIG. 1. The middle six eigenvalues A as a function of 
non-Hermiticity parameter 9 for cf)x = 4>v ~ 0- 

The action corresponding to jl^ ) at 9 — 7r/2 is the 
Wick rotation of the action at = 0. A natural ques- 
tion is whether this analytic continuation to imaginary 
time is accompanied by a smooth variation in properties 
of eigenstates, or whether a phase transition, e.g., metal 
to insulator, occurs as a function of 9 (see ref. [Q). For 
finite 0, we have shown that the geometry of eigenvalue 
degeneracies can change dramatically, but the notion of 
Chern number is still relevant. Whether the entangle- 
ment of eigenvalues can lead to a phase transition in the 
thermodynamic limit is unknown. 

In summary, we extend the notion of Chern number 
to non-Hermitean operators and find that it remains 
a useful tool for distinguishing extended and localized 
states. Entanglement of eigenstates can lead to an intri- 
cate multi-sheeted analytic structure for the Berry flux, 
and the Chern number becomes defined only for the 
group of entangled states as a whole. 
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[7] In the Hermitean case it can be shown that this same 
dynamics is generated by L(X(et)), after subtracting the 
time integral of iJm(X(f)), in the adiabatic limit, e ^ 
but a — et finite. In this limit, tunneling to other eigen- 
states becomes vanishingly small. In the non-Hermitean 
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FIG. 2. Plot of imaginary part of the Berry flux 
lmF±[<j)x, (jiy) for two values of 6, in a small part of the Bril- 
louin zone, [0.357r, 0.407r] x [1.247r, 1.257r]. The Chern numbers 
for eigenvalues I and II are -1-1 and —1 respectively until the 
first branch cut (not shown) occurs at about 6 = .02l7r. As 
shown on the left, another branch cut (the obvious vertical 
lines in the two lower panels separating white from black) 
emerges between 6 = .0257r and 6 = .0297r. The linear grey 
scale spans the range [—1.3, 1.3]. The log-log plot on the right 
displays the divergence of the magnitude of the flux against 
the distance r from an endpoint of the branch cut along the 
paths (h) and (v); the solid line has slope -3/2; deviations 
from it at small r arise from finite numerical resolution. 



case this fails in general because Em is no longer real, 
and an initially infinitesimal amount of tunneling may 
grow exponentially instead of simply oscillating. For our 
purposes the connection with the adiabatic limit is not 
necessary. It matters only that parallel tranport defined 
by (^ along an arbitrary curve C be implemented dy- 
namically by some choice of L(f). 

[8] We assume henceforth that X is a 3-dimensional vector. 

[9] Since 7m (Si) and 7™ (52) differ, Am must have a line 
singularity that emerges from the singularity in Fm and 
passes through Si or S2 ■ This is not a true singularity, but 
a coordinate singularity, as occurs in polar coordinates at 
the poles of a sphere, and its position varies with gauge 
choice. These singularities do not appear in since it 
is gauge invariant. 
[10] Let | — ) is the eigenvector obeying Lmnj — ) = -E|— ). At a 
generic degeneracy one has |-|-) = | — ), and there exists an 
independent generalized eigenvector (also called a cyclic 
vector) \g) that obeys limn\g) ~ | — ) and may be chosen 
so that ((?|±) — 0. For details see an advanced text on 
linear algebra. 

[11] To prove this, observe that the Berry fiux through a 
branch surface is identical for the two eigenvalues in- 
volved, but the surface normals are opposite. Branch 
surfaces make no contribution and the total integral is 
conserved and equal to the sum of the Chern numbers 
before the Chern loop pierced the integration surfaces. 
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